Electromotive force and large-scale magnetic dynamo in a turbulent flow with a mean 

shear 
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An efTect of sheared large-scale motions on a mean electromotive force in a nonrotating turbulent 
flow of a conducting fluid is studied. It is demonstrated that in a homogeneous divergence-free 
turbulent flow the a-effect does not exist, however a mean magnetic field can be generated even in a 
nonrotating turbulence with an imposed mean velocity shear due to a new "shear-current" effect. A 
mean velocity shear results in an anisotropy of turbulent magnetic diffusion. A contribution to the 
electromotive force related with the symmetric parts of the gradient tensor of the mean magnetic 
field (the K-effect) is found in a nonrotating turbulent flows with a mean shear. The K-effect and 
turbulent magnetic diffusion reduce the growth rate of the mean magnetic field. It is shown that a 
mean magnetic field can be generated when the exponent of the energy spectrum of the background 
turbulence (without the mean velocity shear) is less than 2. The "shear-current" effect was studied 
using two different methods: the r-approximation (the Orszag third-order closure procedure) and 
the stochastic calculus (the path integral representation of the solution of the induction equation, 
Feynman-Kac formula and Cameron-Martin-Girsanov theorem). Astrophysical applications of the 
obtained results are discussed. 

PACS numbers: 47.65.-|-a; 47.27.-i 



I. INTRODUCTION 

Generation of magnetic fields by a turbulent flow of a 
conducting fluid is a fundamental problem which has a 
number of applications in solar physics and astrophysics, 
geophysics and planetary physics (see, e.g., [ESHS 
|3]). It is known that small-scale magnetic fluctuations 
with a zero mean magnetic field can be generated in a 
homogeneous nonhelical and nonrotating turbulence by 
a stretch-twist-fold mechanism (see, e.g., |6ll7lla. 19. llOl 
IT]|'). On the other hand, in a homogeneous divergence- 
free turbulent flow the helicity and the a-effect vanish. 

However, the mean magnetic field can be generated in 
a rotating homogeneous turbulent flow due to the com- 
bined action of the fl x J-effect and a nonuniform (dif- 
ferential) rotation [H El El El El, where ft is the 
angular velocity and J is the mean electric current. The 
evolution of the mean magnetic field B is determined by 
equation 



- = Vx(UxB. 



£-DraVxB) 



(1) 



where U is the mean velocity, is the magnetic diffu- 
sion due to the electrical conductivity of fluid, £ = (uxb) 
is the mean electromotive force. The general form of the 
mean electromotive force was suggested in using the 
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symmetry arguments: 

= ayS, -/3y(VxB), + (V(°ff)xB), 

-[^x(VxB)],-Kyfc(9S),fc , (2) 

where {dB)ij — {ViBj + Vji?i)/2, u and b are fluctua- 
tions of the velocity and magnetic field, respectively, an- 
gular brackets denote averaging over an ensemble of tur- 
bulent fluctuations, the tensors a^- and (3ij describe the 
a-effect and turbulent magnetic diffusion, respectively, 
is the effective diamagnetic (or paramagnetic) ve- 
locity, Kijk and 5 describe a nontrivial behavior of the 
mean magnetic field in an anisotropic turbulence. The 
r2 X J-effect, e.g., is associated with the J-term in the 
mean electromotive force. 

In the present paper we suggested a new mechanism 
of generation of a mean magnetic field by a nonrotating 
homogeneous and nonhelical turbulence with an imposed 
mean velocity shear. This mechanism is associated with 
a " shear-current" effect determined by the ^-term in the 
mean electromotive force. On the other hand, the turbu- 
lent magnetic diffusion and the K-effect can reduce the 
growth rate of the mean magnetic field. The K-effect 
arises in an anisotropic turbulence caused by the mean 
velocity shear. Our analysis of the mean-field magnetic 
dynamo showed that the generation of a mean magnetic 
field can occur when q < 2, where q is the exponent 
of the energy spectrum of the background turbulence 
(without a mean velocity shear). In particular, in Kol- 
mogorov background turbulence with g = 5/3 a mean 
magnetic field can be generated. The "shear-current" 
effect was studied using two different methods: the r- 
approximation (the Orszag third-order closure procedure 
[lq| , see Section IV) and the stochastic calculus (the path 



2 



integral representation of the solution of the induction 
equation, Feynman-Kac formula and Cameron-Martin- 
Girsanov theorem, see Appendixes A and B). We also 
calculated the mean electromotive force for an arbitrary 
weakly inhomogeneous turbulence with an imposed mean 
velocity shear. The inhomogeneity of turbulence and 
mean velocity shear cause the a-effect and the effective 
drift velocity of the mean magnetic field. 

The (5-term in the electromotive force which is respon- 
sible for the " shear-current" effect was also independently 
found in j2| in a problem of a screw dynamo using the 
modified second-order correlation approximation. Note 
also that for homogeneous and nonhelical flows another 
mechanism for the magnetic dynamo associated with a 
"negative turbulent magnetic diffusivity" was recently 
discussed in [23,|2ll|23. 

This paper is organized as follows. In Section II the 
general form of the mean electromotive force which in- 
cludes the shear-current effect is obtained using simple 
symmetry reasoning, and the mechanism for the shear- 
current effect is also discussed. In Section III the govern- 
ing equations for turbulent velocity and magnetic fields 
are formulated which then are used for study an effect of 
a mean velocity shear on a turbulence (Section IV) and 
on a cross-helicity (Section V). This allows us to deter- 
mine the mean electromotive force in a turbulent flow of 
a conducting fluid with an imposed mean velocity shear 
(Section VI). The implications of the results for the mean 
electromotive force to the mean-field magnetic dynamo in 
a nonrotating homogeneous turbulence are performed in 
Section VII. Conclusions and astrophysical applications 
of the obtained results are discussed in Section VIII. In 
Appendixes A and B the "shear-current" effect is stud- 
ied using another approach, i.e., stochastic calculus. In 
Appendix C we derived identities used for the derivation 
of equations for the second moment of the velocity field 
and the cross-helicity tensor. 



II. THE QUALITATIVE DESCRIPTION 

The mean electromotive force can be written in the 
form: 



= a,jBj + h,,kBj,k + 0{V^B,) . (3) 

Following to we use an identity Bj,k = {dB)jk — 
£jki{'^'x'B)i/2 which allows us to rewrite Eq. (|2Jl for the 
mean electromotive force in the form of Eq. |(2Jl, where 
Eijk is the fully antisymmetric Levi-Civita tensor, and 

aij = {aij + aji)/2 , (4) 
Vj:'""^ = efc,,ay/2 , (6) 

5^ = {b,j,-h,,j)/A, (7) 
KyTc = -{bijk +bikj)/2 . (8) 



In a homogeneous and nonhelical turbulence the tensor 
Uij vanishes, which implies that aij = and V^'^^'^ = 0. 
Below we consider this case. 

The general form of the mean electromotive force in 
a turbulent flow with a mean velocity shear can be ob- 
tained even from simple symmetry reasoning. Indeed, 
the electromotive force £ is a true vector, whereas the 
mean magnetic field B is a pseudo- vector. Therefore, 
the tensor bijk is a pseudo-tensor (see Eq. (jSJ). For 
homogeneous, isotropic and nonhelical turbulence the 
tensor b^jk = /S^eijk, where /J^ is the turbulent mag- 
netic diffusion coefficient. In a turbulent flow with an 
imposed mean velocity shear, the tensor bijk depends 
on the true tensor \/ jtJi. Note that the tensor VjC7i 
can be written as a sum of the symmetric and anti- 
symmetric parts, i.e., V jUi = {SU)ij - {l/2)£^jkWk, 
where {SU)ij — C^iUj + \7jUi)/2 is the true tensor and 
W = VxU is the mean vorticity (pseudo- vector). Here- 
after we take into account the effect which is linear in 
W jUi. Therefore, the pseudo- tensor bijk has the follow- 
ing general form 

m 

mk + 0,2 £ik {5U)mj 
+03 £jkm{5(l),ni + 04 % I^fc + 05 5ikWj] , (9) 

where ak are the unknown coefficients, Iq is the maxi- 
mum scale of turbulent motions, and the term cx SjkWi 
vanishes since V • B = (see Eq. ||2J)). Using Eqs. 
© we determine the turbulent coefficients defining the 
mean electromotive force for a homogeneous and nonheli- 
cal turbulence: 

= (3^S,j-2(3ollidU),j , (10) 
S = llSoW , (11) 

Kijk = ll (k1 Sij Wk + K2 Sijm idU)mk) , (12) 

where 

Pa = (01-02-203)74, Jo = (04 - 05)72 , (13) 

Kl = -(04 + 05), K2 = -(01+02), (14) 

and (3^ = uqIq/S is the coefficient of isotropic part of 
turbulent magnetic diffusion, while the second term in 
Eq. (|10|l determines anisotropic part of turbulent mag- 
netic diffusion caused by the mean velocity shear. Here 
Mo is the characteristic turbulent velocity in the maxi- 
mum scale of turbulent motions. The K-effect (deter- 
mined by the tensor Kijk) describes a contribution to 
the electromotive force related with the symmetric parts 
of the gradient tensor of the mean magnetic field and 
arises in an anisotropic turbulence caused by the mean 
velocity shear. Since the tensor Kijk is multiplied by the 
symmetric tensor (dB)jk in the the mean electromotive 
force, this allows us to rewrite the tensor Kijk determined 
by Eq. H12I) in a more simple but not in a symmetric 
form. We will show in this paper that the <5-term in 
Eqs. (|2Jl and 111|) for the mean electromotive force de- 
scribes the "shear-current" effect which can cause the 
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mean-field magnetic dynamo in a homogeneous nonro- 
tating turbulence with an imposed mean velocity shear. 

Consider a homogeneous divergence-free turbulence 
with a mean velocity shear, e.g., U — {0,Sx,0) and 
W — (0,0,5"). The mean magnetic field is determined 
by equation: 



dB 

'dt 



Vx[UxB-/3(VxB) 

-6x{VxB)-K{dB)] , 



(15) 



where /3 = Ptj and k = Hijk. Now for simplicity 
we use the mean magnetic field in the form B ~ 
{B^{z),By{z),0). Then Eq. (O reads 



dt 

dBy 

dt 



— S ll (To By + B'J. 
S Bx+ Pj. By , 



where 



ctq = <5o - /3o - Ki/2 - K2/4 



(16) 
(17) 

(18) 



and B" = d'^Bjdz^. In Eq. (flTjl we took into account that 
Z§_B" ^ B^^ i.e., the characteristic spatial scale Lb of the 
mean magnetic field variations is much larger than the 
maximum scale of turbulent motions Iq. This assumption 
corresponds to the mean-field approach. The first term 
(oc SBx) in RHS of Eq. jT7|l_ r)lavs a role of the differential 
rotation. Indeed, Vx(UxB) = (B • V)U - (U • V)B = 
SBx^y, and for the chosen configuration of the mean 
magnetic field, (U • V)B = 0. 

A solution of Eqs. (|16|l and H17|) we seek for in the form 
oc exp(7f -|- iKz), where 7 is given by 



(19) 



where ctq = (02 -t- 03 -|- 2ai)/2. The magnetic dynamo 
instability can be excited when ctq > 0. In this paper we 
will find unknown coefficients a^, which will allow us to 
determine the conditions for the generation of the mean 
magnetic field due to the magnetic dynamo instability 
caused by the "shear-current" effect. 

In order to elucidate the physics of the " shear-current" 
effect, let us compare the a-term in the electromotive 
force which is responsible for the generation of the mean 
magnetic field, i.e.. 



(20) 



(see, e.g., 0,^3)' '^i^h the 5-term in the electromotive 
force caused by the "shear-current" effect, i.e., 



£f = -(^x(VxB)), oc -(W • V)B, 



(21) 



where A = V(u^)/ (u^) determines the inhomogeneity of 
turbulence. Here for simplicity we considered an isotropic 
Qf-tensor, i.e., a^- = a6ij. There is an analogy between 
the a-term and the 5-term in the electromotive force. In 



particular, the mean vorticity W plays a role of rota- 
tion rj and an inhomogeneity of the mean magnetic field 
plays a role of the inhomogeneity of turbulence in the 
afi-dynamo (see below). During the generation of the 
mean magnetic field in both cases, the mean electric cur- 
rent along the original mean magnetic field arises. The a- 
effect is related with the hydrodynamic helicity oc (f2 • A) 
in an inhomogeneous turbulence. The deformation of 
the magnetic field lines is caused by upward and down- 
ward rotating turbulent eddies. Since the turbulence is 
inhomogeneous (which breaks a symmetry between the 
upward and downward eddies), their total effect on the 
mean magnetic field does not vanish and it creates the 
mean electric current along the original mean magnetic 
field. 

In a turbulent flow with an imposed mean velocity 
shear, the inhomogeneity of the original mean magnetic 
field breaks a symmetry between the infiuence of up- 
ward and downward turbulent eddies on the mean mag- 
netic field. The deformation of the magnetic field lines is 
caused by upward and downward turbulent eddies which 
causes the mean electric current along the mean magnetic 
field and produces the magnetic dynamo. The magnetic 
dynamo instability due to the "shear-current" effect is 
determined by a system of Eqs. ifTCI) and ((T7|l . and there 
is a coupling between the components of the mean mag- 
netic field. In particular, the field By generates the field 
Bx due to the "shear-current" effect (see Eq. (dJ). This 
is similar to the a effect. On the other hand, the field B^ 
generates the field By due to the pure shear effect (see 
Eq. HITI) ). like the differential rotation in ail-dynamo. 

In the next sections we will describe the above mag- 
netic dynamo effect quantitatively using two different 
methods: the r-a ppr oximation (the Orszag third-order 
closure procedure [13) and the stochastic calculus (the 
path integral representation of the solution of the in- 
duction equation, Feynman-Kac formula and Cameron- 
Martin-Girsanov theorem). 



III. THE GOVERNING EQUATIONS 

Our goal is to study an effect of sheared large-scale 
motions on a mean electromotive force in a nonrotating 
turbulent flows of a conducting fluid. The momentum 
equation for the fluid velocity v and the induction equa- 
tion for the magnetic field h read 



d \ VP 

-+vV V = + F,^(h) + ^Av + F(''*), 

dt J p 



dt 



V • V h = (h • V)v + D,nAh , 



(22) 
(23) 



where V • v = and V • h = 0, v is the kinematic 
viscosity, Fm(h) = — (l//ip)[hx(Vxh)] is the magnetic 
force, n is the magnetic permeability of the fluid, F^***' 
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is the random external stirring force, P and p are the 
pressure and density of fluid, respectively. 

We will use a mean field approach whereby the veloc- 
ity, pressure and magnetic field are separated into the 
mean and fluctuating parts: v = U + u, P = P + p, and 
h = B + b, the fluctuating parts have zero mean values, 
and U = (v), P = (P), B = (h). Averaging Eqs. ^ 
and H23|l over an ensemble of fluctuations we obtain the 
mean- field equations. In particular, the evolution of the 
mean magnetic field B is determined by Eq. (Q, where 
£ = (u X b) is the mean electromotive force. To de- 
termine the mean electromotive force we use equations 
for fluctuations u(t,r) and h{t,r) which are obtained by 
subtracting equations for the mean fields from the corre- 
sponding equations H22(l and (|23|l for the total fields: 



du 

'dt 

dh 

'dt 

where 



-(tJ • V)u - (u • V)U - ^ + F™(b, B) 

P 



-(U • V)b + (b • V)U - (u • V)B 
+ (B- V)u-KBjv , 



(24) 
(25) 



F™(b,B) 

Un 



lip 



[bx(VxB) + Bx(Vxb)] , 



((u • V)u) - (u • V)u + [(bx (Vxb)) 
-bx(Vxb)]/(^p) + j/Au , 
Vx(uxb - (uxb)) + AnAb . 



We consider a turbulent flow with large hydrodynamic 
(Re = IoUq/u ^ 1) and magnetic (Rm = IqUq/D^ ^ 1) 
Reynolds numbers, where Uq is the characteristic veloc- 
ity in the maximum scale Iq of turbulent motions. In 
the next sections we will use Eqs. (|^ and (^5)1 to study 
an effect of a mean velocity shear on a turbulence (Sec- 
tion IV) and on a cross-helicity (Section V) in order to 
determine the mean electromotive force. 



IV. EFFECT OF A MEAN VELOCITY SHEAR 
ON A TURBULENCE 

In this section we study an effect of a mean velocity 
shear on a turbulence using Eq. (|24|) . We neglect an 
effect of the mean magnetic field on the turbulence, i.e., 
we neglect the magnetic force Fm(b, B) in Eq. ll^ . This 
is valid when / p, <C p{u'^)/2, i.e., we do not consider 
the quenching effects (see, e.g., [H HI IH IH 18] ) . 
We use a two-scale approach, i.e., a correlation function 
is written as follows 



(ii,(x)uj(y)) 



Wi(ki)iij(k2)) exp[i(ki-x 
k2-y)] dki dk2 
/ij (k,R) exp(ik-r) dk , 



/y(k,R) = J{u,Ck + K/2)u,{-'k- 



-K/2)) 



X exp(iK-R) dK 



(see, e.g., |23, ISfl), where R and K correspond to the 
large scales, and r and k to the small scales, i.e., R — 
(x-Ky)/2, r = x-y, K = ki ka, k = (ki - 
k2)/2. We assume that there exists a separation of scales, 
i.e., the maximum scale of turbulent motions Iq is much 
smaller then the characteristic scales of inhomogeneities 
of the mean fields. 
Now we calculate 



9/,j(ki,k2) 
dt 



= (-Pm(ki 



,(k2)) 



(u,(ki)P,„(k2) 



dt "^^ 

5u„(k2) , 



dt 



(26) 



where we multiplied equation of motion H24|l rewritten 
in k-space by Pij (k) — Sij — fcy in order to exclude the 
pressure term from the equation of motion, 6ij is the 
Kronecker tensor and fcy — kikj/k'^. Thus, the equations 
for /ij (k, R) is given by 



9fv 
dt 



tjmn 



(N) 



(27) 



where 



'-ijmn 



(U) 



2 q Sim Spn ^iT7i Sj q Sjip 



SimSjnkq-g^ ) VpJ7, , (28) 



and 
ing due 



(N) 



(k, R) is the third moment appear- 
to the nonlinear term, V = d/dR, 



F,,(k,R) = (J^,(k,R)u,(-k,R)) + (u,(k,R)F,(-k,R)) 
andF(k,R,i) = -kx (kxF("*)(k, R))/fc2. Equation j23 
is written in a frame moving with a local velocity U of 
the mean flows. In Eqs. H27|) and (|28|l we neglected small 
terms which are of the order of 0(|V^U|). Note that 
Eqs. (|27l) and H28|) do not contain terms proportional to 
OdV^tJI). To get Eqs. ^ and ^ we used an identity 
derived in Appendix C. 

Equation H27() for the second moment (k, R) con- 
tains third moment f^j^'^ (k, R) and a problem of closing 
the equations for the higher moments arises. Various ap- 
proximate methods have been proposed for the solution 
of problems of this type (see, e.g., |3 |^ H^l)- The 
simplest procedure is the r-approximation (the Orszag 
third-order closure procedure [ill)- For magnetohydro- 
dynamic turbulence this approximation was used in [s^ 
(see also |2^ |3 113 )■ In the simplest variant, it al- 
lows us to express the deviations of the third moment 



moment fij (k, R 



/j " (k, R) in terms of that for the second 



(AT) 



/i°^(k,R): 



ANO) 



fij fq 



(0) 



(29) 
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where the superscript (0) corresponds to the background 
turbulence (it is a turbulence with zero gradients of the 
mean fluid velocity, ViC/j = 0), and r(/c) is the corre- 
lation time of the turbulent velocity field. Here we as- 
sumed that the time T{k) is independent of gradients of 
the mean fluid velocity because in the framework of the 
mean- field approach we may only consider a weak shear: 
To|VJ7| < 1, where tq = Iq/uq. 

The T-approximation is in general similar to Eddy 
Damped Quasi Normal Markovian (EDQNM) approxi- 
mation. However some principle difference exists between 
these two approaches (see [li 0). The EDQNM clo- 
sures do not relax to equilibrium, and this procedure does 
not describe properly the motions in the equilibrium state 
in contrast to the T-approximation. Within the EDQNM 
theory, there is no dynamically determined relaxation 
time, and no slightly perturbed steady state can be ap- 
proached In the T-approximation, the relaxation 
time for small departures from equilibrium is determined 
by the random motions in the equilibrium state, but 
not by the departure from equilibrium |18| |. As follows 
from the analysis by |l8l | the T-approximation describes 
the relaxation to equilibrium state (the background tur- 
bulence) much more accurately than the EDQNM ap- 
proach. 

Note that we applied the T-approximation 129() only 
to study the deviations from the background turbulence 
which are caused by the spatial derivatives of the mean 
velocity. The background turbulence is assumed to be 
known. Here we use the following model of the back- 
ground nonhelical, isotropic and weakly inhoniogeneous 
turbulence: 



/^^(k,R) = 



1 



^y(k) 



fc,V,) )ulE{k,-R) , 



(30) 



where T(fc) = 2Tof(fc), E{k) = -df{k)/dk, f{k) = 
(fc/fco)^"*, 1 < g < 3 is the exponent of the kinetic en- 
ergy spectrum (e.g., 9 = 5/3 for Kolmogorov spectrum), 
fco = 1/^0- 

We assume that the characteristic time of variation of 
the second moment fij (k, R) is substantially larger than 
the correlation time T{k) for all turbulence scales. Thus 
in a steady-state Eq. (|27|l reads 



[5^„^5,,^ - t4„,„ (U)] [/„„ (k, R) - / ^ ] 



= tL. 



(0) 



(31) 



where we used Eq. (j^HJ- The term Fij in Eq. (|77|l de- 
termines the background turbulence. The solution of 
Eq. H31|) yields the second moment fij (k, R) : 

(k, R) « /(°) + t4„„ (U)/i«) (k, R) , (32) 

where we neglected terms which are of the order of 
0(|tVUP) <C 1. The first term in Eq. ^ is indepen- 



dent of the mean velocity shear and it describes the back- 
ground turbulence. The second term in Eq. (|32|l deter- 
mines an effect of the mean velocity shear on turbulence. 



V. EFFECT OF A MEAN VELOCITY SHEAR 
ON THE CROSS-HELICITY 

In order to determine the mean electromotive force 
£i(r = 0, R) = Einm J ffmn(k, R)dk, we derive equation 
for the cross-helicity tensor: 

g„(k,R) = y"(&,(k + K/2)7i,(-k + K/2)) 
X exp(iK-R) dK , 
using Eqs. 1)241) and (|25|l . i.e., we calculate 

%y(ki,k2) _ 5u„(k2) 

= {bi[ki)Pjn[K2)- ^ 



dt 



at 



+ ( %(k2)) 



(33) 



This yields equation for (7y (k, R), 

— '^ijmn{'^)9mn ~^ ^ijmni^) fmn ~t~ 9ij : (34) 



which describes an effect of a mean velocity shear on the 
cross-helicity, where 

t^zjmn(U) — Ik j q5 i^fyi^ fiU q ^itn^nUj -|- ^jn^mUi 

d 

+<5im(5j„(Vg[/p)fcp^^ , (35) 



-^^zjmn(-B) — ^im^jn-Bp ^^^p ~t~ ) ^jn-Bi 

2 ^ira^ jn^p.q'^p 



(36) 



and g\^^ (k, R) is the third moment appearing due to 
the nonlinear terms, Bi,j = dBi/dRj. Equation H34I) is 
written in a frame moving with a local velocity U of the 
mean flows. To get Eqs. (jS^-JSHl) we used an identity 
derived in Appendix C. 

Now we use the T-approxiniation which allows us to 
express the third moments f/|j^'*(k, R) in terms of the 
second moments gij (k, R) : 



(N) 

9ij 



9ij 

m 



(37) 



where f (fc) is the characteristic relaxation time, and we 
took into account that the cross-helicity tensor gij for 
B = is zero, i.e., gijCB = 0) = 0. Note that we apphed 
the T-approximation l|37() only to study the deviations 
from the original turbulence (i.e. the turbulence with 
B = 0). These deviations are caused by a weak mean 
magnetic field. We considered the case when the original 
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turbulence does not have magnetic fluctuations. Now we 
assume that the characteristic time of variation of the 
second moment gij (k, R) is substantially larger than the 
correlation time T{k) ^ T{k) for all turbulence scales. 
This allows us to get an equation for the cross-helicity 
tensor gtj (k, R) in a steady state: 

5?7in (k, R) [Jirnt^jn '^<^ijmn(U)] = TAIijmni^'^^ fnin ;(38) 

where /m„(k, R) in Eq. (jSHl is determined by Eq. (p!^ . 
The solution of Eq. (jSHJ yields 

5y (k, R) = TMy „„ (B)/i°^ + t2 [Mijmn (B)/„„Ufc (U) 



+ Jijmn{'^)Mmnlk(B)]fl 



(0) 
Ik ' 



(39) 



where we neglected terms which are of the order of 
OdrVUp) < 1. The first term in Eq. ^ is inde- 
pendent of the mean velocity shear and it describes the 
isotropic turbulent magnetic diffusion. The second term 
in Eq. (|39|l determines an effect of the mean velocity shear 
on turbulence which causes a modification of the mean 
electromotive force, i.e., this term describes an indirect 
modification of the mean electromotive force. The last 
term in Eq. H39|l determines a direct modification of the 
mean electromotive force by the mean velocity shear. 



VI. THE MEAN ELECTROMOTIVE FORCE 

Using Eqs. H32|) and (|39|) we determine the mean elec- 
tromotive force for an inhomogeneous turbulence with a 
mean velocity shear: 

£i(r = 0,R) = Einrn J g„in{K'R) dk 

= UijBj +bi.jkBj^k , (40) 



where 



/2 
18 



6A, 



Aq-9 



+3e^Jm{WxAU + (5y(W • A) - mAjj , (41) 

f3^£r,k - ^[2(8g - ll)e^j„ddU)^k 

+ 18eifcm(9C/)„j + (4g - n)5,jWk 

+ {Aq-7)5^kWj]. (42) 

Here A = V(u^)/(u^) and /q — t^uq. Equations lUT)) and 
(|42|l allows us to obtain the turbulent coefficients defining 
the mean electromotive force: 

/2 / 1 
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-JL[ (W.A)%--(W^,A,+W^,A,, 



[dU] {dUU 

(43) 

V(°^) = |('(49 + l)(aC/).,A, + |(WxA)') , (44) 



and the tensor of turbulent magnetic diffusion , the 
^-effect and the tensor k^j; are determined by Eqs. H1(J|) - 
(O with 

/3o - 2(5 - 2g)/45 , <5o = 1/9 , (45) 
= -8(3-g)/45, K2 = 4(49- l)/45 , (46) 

where we used Eqs. Q-ljHIl. It is seen from Eqs. H43|l 
and (|44|) that the a-effect described by the tensor ay 
and the effective drift velocity V(°^^ of the mean mag- 
netic field require an inhomogeneity of turbulence (i.e., 
A 7^ 0). The K-effect determined by the tensor K^fc arises 
in an anisotropic turbulence caused by the mean veloc- 
ity shear. We will show in the next section that the 
5-term in the equation for the mean electromotive force 
can cause the mean-field magnetic dynamo in a homo- 
geneous nonrotating turbulence with an imposed mean 
velocity shear. 



VII. THE MEAN-FIELD MAGNETIC DYNAMO 
IN A HOMOGENEOUS TURBULENCE WITH A 
MEAN SHEAR 

Consider a homogeneous divergence-free turbulence 
with a mean velocity shear, e.g., tj = (0,53;, 0). In this 
case W = (0,0,5), A = 0, the a-effect and the effec- 
tive drift velocity V^®'^) of the mean magnetic field van- 
ish. The mean magnetic field is determined by Eq. H15|l . 
where 5 = (Jo^qW describes the "shear-current" effect 
and j3ij = P^Sij — 2(3olQ{dU)ij corresponds to the turbu- 
lent magnetic diffusion with an anisotropic part oc /3o- 
The tensor Kijk = ll{Ki5ijWk + K2£ijm idU)mk) de- 
scribes a contribution to the electromotive force related 
with the symmetric parts of the gradient tensor of the 
mean magnetic field and arises in an anisotropic turbu- 
lence caused by the mean velocity shear. Since the tensor 
Kijk is multiplied by the symmetric tensor {dB)jk in the 
the mean electromotive force, this allows us to rewrite 
the tensor Kijk in a more simple but not in a symmetric 
form. For simplicity we use the mean magnetic field in 
the form B = {Bj:{z), By{z),0). Then Eq. ((THJl reduces 
to Eqs. (|15|l and (|17|l . where the parameters Sq, (3o, ki 
and K2 are determined by Eqs. H45|l and (|46|l . 

A solution of Eqs. 116|) and (|17|) we seek for in the 
form cx eicp{'jt + iKz). Thus the growth rate of the mean 
magnetic field is given by 



J === SIqK y/6a- i3q - kq - K"^ 



(47) 



where kq = (2ki -|- K2)/4 (817 - 13)/45. It follows 
from Eq. H47() that the "shear-current" effect (cx Sq) 
causes the generation of the mean magnetic field, whereas 
the anisotropic (cx /3o) and isotropic (cx [3^) turbu- 
lent magnetic diffusions and the K-effect (cx kq) re- 
duce the growth rate of the mean magnetic field. Note 
that the maximum growth rate of the mean magnetic 
field, 7inax — S'^lo{5{) ~ Po — Ko)/4/3y, is attained at 
K = Km = SIq^/Sq — Po — '«o/2/3t. Using expressions 
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for 5q , Pq and kq we rewrite the growth rate of the mean 
magnetic field in the form 

l^K^^^Y'sioK-P^K^. (48) 

Therefore, the generation of the mean magnetic field is 
possible when the exponent q of the energy spectrum 
of the background homogeneous turbulence (without im- 
posed mean velocity shear) is less than 2. Thus, in the 
Kolmogorov background turbulence with q = 5/3 the 
mean magnetic field can be generated due to the " shear- 
current" effect. The sufficient condition 7 > for the dy- 
namo instability reads Lb/Iq > ny/5/{ToS\/2 — q), where 
Lb = 2tt/K. 

The magnetic dynamo instability due to the "shear- 
current" effect is different from the magnetic instability 
suggested in The latter instability is caused by the 
"negative turbulent magnetic diffusivity" and is deter- 
mined by Eq. (4.2) for in This equation is de- 
coupled from that for the field By, i.e., there is no real 
coupling between the components of the mean magnetic 
field Bx and By. In contrast to this, the magnetic dy- 
namo instability due to the "shear-current" effect is de- 
termined by a system of equations (|16|l and (|17() for the 
components B^ and By. This implies that there is a cou- 
pling between these components of the mean magnetic 
field. In particular, the field By generates the field B^ 
due to the (5-term ("shear-current" effect), see the first 
term in Eq. H16|l. This is similar to the a effect. On 
the other hand, the field Bx generates the field By due 
to the pure shear effect (see the first term in Eq. (EJ), 
like the differential rotation in all-dynamo. In this sense 
the instability due to the "shear-current" effect is a pure 
magnetic dynamo instability. 

However, the above mechanism is different from that 
for afi-dynamo. Indeed, the dynamo mechanism due 
to the "shear-current" effect acts even in homogeneous 
small-scale turbulence, while the alpha effect vanishes for 
homogeneous turbulence. The difference between these 
magnetic dynamo mechanisms can be seen in the form 
of the growth rate of the mean magnetic field. Indeed, 
the generation of the mean magnetic field is caused by 
a coupling between the "shear-current" effect (described 
by the first term in Eq. which is proportional to 

the second-order spatial derivative of the mean magnetic 
field) and the pure shear effect (described by the first 
term in Eq. H17|l which is proportional to the mean mag- 
netic field). Then the first term in the expression for 
growth rate in Eq. (|47|) (which is responsible for the gen- 
eration of the mean magnetic field due to the "shear- 
current" effect) is proportional to the wave number K . 

On the other hand, the afi-dynamo is caused by a 
coupling of the a-effect (the corresponding term in the 
mean-field equation is proportional to the first-order spa- 
tial derivative of the mean magnetic field) and the differ- 
ential rotation (the corresponding term is proportional 
to the mean magnetic field). Then the term in the ex- 
pression for growth rate of the instability (which is re- 



sponsible for the generation of the mean magnetic field 
in the ai7-dynamo) is proportional to K^^^. 

Note that the properties of the magnetic dynamo 
caused by the "shear-current" effect are also different 
from that for the rjxj-effect. In particular, the mean 
magnetic field can be generated due to the JlxJ-effect 
for an arbitrary exponent q of the energy spectrum of 
the background homogeneous turbulence (see dfi]). The 
Ox J~effect is caused by the term 6x (VxB) in the mean 
electromotive force. In a slow rotating {flro <C 1) and 
homogeneous turbulence S = —{2/9)11^1 (for details, see 
[3|)- Note also that the Ox J-effect cannot generate the 
mean magnetic field without a differential rotation. 

VIII. CONCLUSIONS 

In the present paper we discussed a new mechanism 
of a generation of a mean magnetic field by a nonro- 
tationg and nonhelical homogeneous turbulence with an 
imposed mean velocity shear. This mechanism is asso- 
ciated with a "shear-current" effect. We showed that 
when the exponent of the energy spectrum of the back- 
ground turbulence (without the mean velocity shear) is 
less than 2, a mean magnetic field can be generated. 
We calculated the mean electromotive force for an arbi- 
trary weakly inhomogeneous turbulence (A/q <C 1) with 
an imposed mean velocity shear. Inhomogeneity of tur- 
bulence and mean velocity shear cause the a-effect and 
the effective drift velocity of the mean magnetic field. 
The " shear-current" effect was studied using two different 
methods: the r-approximation (the Orszag third-order 
closure procedure) and the stochastic calculus (the path 
integral representation of the solution of the induction 
equation, Feynman-Kac formula and Cameron-Martin- 
Girsanov theorem, see Appendixes A and B). 

The obtained results may be important in astrophysics, 
e.g., in extragalactic clusters and in interstellar clouds. 
The extragalactic clusters are nonrotating objects which 
have a homogeneous turbulence in the center of a ex- 
tragalactic cluster. Sheared motions are created between 
interacting clusters. The observed magnetic fields cannot 
be explained by a small-scale turbulent magnetic dynamo 
(see, e.g., @). It is plausible to suggest that the "shear- 
current" effect can produce a mean magnetic field in the 
extragalactic clusters. The sheared motions can be also 
formed between interacting interstellar clouds. The lat- 
ter can result in a generation of a mean magnetic field. 
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APPENDIX A: INVESTIGATION OF THE 
"SHEAR-CURRENT" EFFECT USING 
STOCHASTIC CALCULUS 

In this Appendix we study the "shear-current" effect 
using stochastic calculus for a random velocity field with 
a finite correlation time. In order to derive an equation 
for the mean magnetic field we use an exact solution of 
the induction equation H23|l in the form of a functional 
integral for an arbitrary velocity field taking into account 
a small yet finite magnetic diffusion caused by the elec- 
trical conductivity of fluid. This magnetic diffusion, Dm, 
can be described by a random Brownian motions of a par- 
ticle. The functional integral implies an averaging over 
a random Brownian motions of a particle. The form of 
the exact solution used in the present paper allows us 
to separate the averaging over both, a random Brownian 
motions of a particle and a random velocity field. This 
method yields the solution of the induction equation H23|) 
with an initial condition h(< = s, x) = h(s, x) in the form 

/i,(t,x) = M^{G„(t,5,0 exp(r • V)/ij(s,x)} (Al) 

(see Appendix B and 0), where ^* = ^ — x, 
Gy (i, s,^) is determined by equation dGij{t, s,$^) / ds = 
NikGkj (t, s, with the initial condition Gij {t — s) — Stj . 
Here Nij = dvi/dxj, M^{-} denotes the mathematical 

expectation over the Wiener paths ^ = x — ^ v(i — 
a, ^) da + (2Z)m)^/^vif(i — s), and the magnetic diffusion, 
Djn, is described by a Wiener process w{t). 

Consider a random velocity field with a finite con- 
stant renewal time. Assume that in the intervals 
. . . (— r, 0]; (0, r]; (t, 2t]; . . . the velocity fields are statis- 
tically independent and have the same statistics. This 
implies that the velocity field looses memory at the pre- 
scribed instants t = kr, where fc = 0, ±1, ±2, .... This 
velocity field cannot be considered as a stationary veloc- 
ity field for small times ~ r, however, it behaves like a 
stationary field for t ^ t. 

In Eq. ljAl|l we specify instants t = {m + l)r and 
s — rriT. Note that the fields hj(mT,x.) and Gij{{m + 
l)r, mr, ^) are statistically independent because the field 
hj{mT,:x.) is determined in the time interval (— oo,TOr], 
whereas the function Gij{{jn + I)t, mr, ^) is defined 
on the interval {mT,{m -\- 1)t]. Due to a renewal, the 
velocity field as well as its functionals hj {mT, x) and 
Gy ((m-|- 1)t, mr, ^) in these two time intervals are statis- 
tically independent. Averaging Eq. ljAl|l over the random 
velocity field yields the equation for the mean magnetic 
field 

S,((m+l)T,x) = M^{(G,,(t,s,0 exp(r • V))} 

xi?j(mT,x), (A2) 
where the operator exp(^* • V) is determined by 

1 



exp(r -V) = l + C -V + -Vf 

+^(r •vr + ... , 

to! 



(A3) 



and the angular brackets (•) denote the ensemble av- 
erage over the random velocity field. Note that |(^* • 
V)B|/|B| Iq/Lb < 1. Thus in the framework of the 
mean-field approach we can neglect in Eqs. (|A2|) and (|A3p 
the terms ^ 0[(^* • V)'^B]. Now we use the identity 



B,{t + T,yL) = exp^T^jB,(i,x) , 
which follows from the Taylor expansion 



(A4) 



d 



f^t + r) = Y^[T-] fit)=c^p(r 



d_\m 

dt to! 



Therefore, Eqs 

d_ 
dt 



expl T 



2ll-(IA4ll yield 

Bi{t,X.) — {Gij + GijS^m^r. 



-Gy>„„V™V„)i?j = exp(rL)B ,(A5) 



where Gy = M|{(G„)} = % + C/,,,t + 0[{VU)\ 

^^ = M^{(e*)} = -C7,T + 0[{VU)% C7,,, = dU,/dx„ 

A,,n = ^f|{(G,yC„)}, G.,„„ = M|{(G.yCCJ}, and 

we introduced the operator L, which allows us to reduce 
the integral equation ljA2|l to a partial differential equa- 
tion. Indeed, Eq. IjASp . which is rewritten in the form 



exp 



T L 



d_ 

dt 



B = B 



reduces to 



9B 

'dt 



= LB 



(A6) 



(A7) 



The Taylor expansion of the function exp(TL) reads 

exp{TL) = E + tL + {tL)^/2 + ... , (A8) 

where E is the unit operator. Thus, Eqs. HA5|I and (|A8() 

yield 



Li-j — -(G., 



D 



Gij^rn^ r 



(A9) 



where £).ym„ = {Gijmn - AikmAkjn) /'^t. Now we con- 
sider homogeneous and nonhelical background turbu- 
lence, then Aijk = and the equation for the mean mag- 
netic field is given by 



dt 



= [V X (UXB)], + AjmnV„V„Bj , (AlO) 



where 



A,>n„ - 2^M^{(G.,CC>} • (All) 

For a turbulent flow with an imposed mean velocity 
gradient, the turbulence is anisotropic. Let us determine 
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the tensor Dijmn in this case. Solution of the equation 
dGij/ds = Ni^Gkj with the initial condition Gy (i = 
s) = Sij is given by 



Jo 



+ f N,k{ts,^)ds f Nkj{ta,^)da + ... , (A12) 
Jo Jo 

which was solved by iterations, where t^- = t + t — a. 
Since the velocity field is separated into the mean and 
fluctuating parts: v = U + u, the tensor Gy can be 
presented in the form 



Gij — gtj + G[j^ 



9tj = 



dxk 



(A13) 

gkjita,^)da, (A14) 



G[f \lJ) = tJ^^k I 9kAt.,i)da 



3 dxk 
Using Eqs. HA12P - HA15I) we obtain 



Gtj ^ gij + gikUk,p / gpj{tcr,$,) da + 0[{VU) ] . 
Jo 

(A16) 

For the derivation of Eq. (jA16p we used the identity: 

(E - X)-^ = E + X + XX + XXX + ...+ , (A17) 

where X is an arbitrary operator. Similarly, the trajec- 
tory ^* can be written in the form: 

C = i^ - U,.k I Ikit.A) da + 0[{VUf] , (A18) 
Jo 

where |i = - Cuiita,^) da + {2Dm)^/'^Wi{T). Using 
Eqs. (jA16p and l|Al8|) we determine the tensor Dij^n 



D, 



where F, 



tjmn 



Uk,p{giklm£,n Jo 9pjda) and K, 



(A19) 



ijmn 



{gij£,m Jq £,n da) . In this section hereafter the angular 
brackets denote the both averaging: the averaging over a 
random velocity field and the averaging over the Wiener 
trajectories. Now we determine the tensor {gij£.m£.n), 
which describes turbulent magnetic diffusion. We take 
into account that in a homogeneous and nonhelical tur- 
bulence without mean shear (Uij = 0), the mean- field 
equation reads 



= -[Vx(/3,VxB)], 



In this case bijk = P^Sijk- Therefore, 

^imk^knj 

= '2Tf3^{SijSmn - StnSmj) , (A21) 



where 



1 r°° 



and we used an identity 



( / u,{^i,i)d^. / Uj{a,i)da) 
Jo Jo 

= 2r/ {u,{QX)uj[a,^))da . (A23) 
The integration of Eq. (|A21|) over r yields the tensor 



(A20) 



Now we construct the tensor Fijmn ■ The general form of 
this tensor reads 

mn mn ) , (A25) 

where we took into account that the tensor Fijmn in 
Eq. I|A10|I is multiplied by a tensor Wm'^nBj which is 
symmetric with respect to indexes (m, n). Since V • B = 
0, the tensor Fijmn does not contain the terms with the 
tensors 6jm and Sjn, and Fijmn satisfies to the condi- 
tion Fijmn^ nB j — 0. The latter equation yields 
Ci = — C3 and C2 — —C4. Equation HA25P does not 
have the term cx SinkSjpm{dU)pk because we considered 
a nonhelical turbulence. This is the reason that Eqs. 142|) 
and (|A27|I does not contain the term oc ejkm{dU)mi (cf. 
with the general form of the tensor bijk determined by 
Eq. ©). Therefore, Eqs. l|IT9)l . l|MT|l . l|X24)l and (|X25|l 
yield the tensor Dijmn ■ 

F^ijmn — /^7^ f^mn '^\G\{Um^j^in UijSmn) 

)-(aC/)„„(5,j]}.(A26) 

Now we use the identities: 

{Um.j^in Ui_j 5mn)VmVnBj = [V X (V X (B • V)U)], 
(dU)mn^rn^nB^ = (V X Q), , 

where Qi = £ijp{dU)pnS/ nBj . The latter identity is de- 
rived using the following identity: £ims£sjp(dU)pn — 
{dU)mnSij — {dU)inSmj- Notc that we neglect the second 
and higher order spatial derivatives of the mean velocity. 
We also neglected the cross-effect terms which describe 
an interaction between molecular and turbulent effects. 
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Thus, Eq. ljA26|) and the above identities allows to deter- 
mine the tensor bijk' 

bijk = f^^eijk + (3^t[{Ci + C2)eikmidU)mj 

+ \{C2 - Ci)%VKfe - e,,,n{diJ)mk] . (A27) 

Using Eqs. Q)-® and HA27|I we determine the turbulent 
coefficients defining the mean electromotive force. They 
are given by Eqs. (CHll-OSl with 

/3o = -^(Ci+C2 + l), ^o = ^(C2-Ci), (A28) 

«i = -i(C2-C'i), n2 = -\{Ci+C2-l), (A29) 

A 



where Iq = ^J^f3~T. For simplicity we use the mean 
magnetic field in the form B = (^^^(z), 0). Then 

Eq. l(T^ reduces to Eqs. and ((T7|l . where the param- 
eters (5o, /3oi and K2 are determined by Eqs. (|A28p 
and HA29() . A solution of Eqs. (fTC|l and lfT7|l we seek for 
in the form cx exp(7t -I- iKz). Thus the growth rate of 
the mean magnetic field due to the shear-current effect 
is given by 



(A30) 



where we used that ctq = ^2/3. Therefore, the magnetic 
dynamo instability can be excited when C2 > 0. 

This approach does not allow us to take into account 
the effect of mean velocity shear on turbulence. The 
method used in this Appendix only describes the effect of 
shear on the cross-helicity tensor. This is one of the rea- 
sons that the results obtained by this method are quan- 
titatively different from that obtained by the r approx- 
imation. However, the form of the electromotive force 
and a possibility for the large-scale magnetic dynamo in 
a homogeneous turbulence due to the shear-current effect 
are clearly demonstrated by two different approaches. 



APPENDIX B: DERIVATION OF EQ. iXTl 

In order to derive Eq. ljAl|l we use an exact solution of 
Eq. (|23|) with an initial condition h{t = s,x) = h(s,x) 
in the form of the Feynman-Kac formula: 

h,{t,^) = M^{G,j{t,s,^{t,s))h,{s,^{t,.s))} , (Bl) 

where the Wiener paths $,{t,s) = x — J* \'[t — 
a, ^(t, cr)] da + (2£'„i)^/'^w(t — s). Now we assume that 

h^{t, €) = y" exp(«C • <l)h,is, q) dq . (B2) 

Substituting Eq. into Eq. jBlll we obtain 

h,{s,x) = y Af|{G,y (i,s,^(t,s)) cxp[i^* • q] /ij (s,q)} 



In Eq. ljB3|l we expand the function exp[i^* • q] in Taylor 
series at q = 0, i.e., exp[i^* • q] = X]fc!Lo(l/'^')(*^* ' Q)*^- 
Using the identity (iq)*^ exp[«x • q] = V*^ exp[ix • q] and 
Eq. |jB3l) we get 



k=0 

X J /ij(s, q) exp(iq • x) dq} . (B4) 



After the inverse Fourier transformation in Eq. ljB4p we 
obtain Eq. (|A1|) . Equation (jB2p can be formally consid- 
ered as an inverse Fourier transformation of the function 
hi{t,^). However, ^ is the Wiener path which is not a 
usual spatial variable. Therefore, it is desirable to derive 
Eq. by a more rigorous method as it is done below. 

To this end we use an exact solution of the Cauchy 
problem for Eq. 1)23(1 with an initial condition h{t = 
s,x) = h(s,x) in the form 

/i,(i,x) - Af^{J(t,s,OG„(t,s,0/ij(s,C(t,s))} , (B5) 

where the matrix is determined by the equation 
dGij{t, s, (^)/ds = NikGkj{t, s, with the initial condi- 
tion Gij{t ~ s) = 6ij, and the function J(t, s, C) is given 
by 



J(i,s,C) = cxp 



1 

V^Dm Jq 



Vi{t - 1],Cit,v))dw^{r]) 



1 



AD. 



m Jo 



w'{t-,^,Cit,ri))dv , (B6) 



X exp(iq • x) dq 



(B3) 



w(i) is a Wiener process, and Af^{-} denotes the math- 
ematical expectation over the paths C(^7 s) — :x. + 
(2D„)i/2(w(t)-w(s)). The solution (|B5l was first found 
in 36]. The first integral * v(i — yy, ^(i, r;)) • d'w{ri) in 
Eq. (|B6)I is the Ito stochastic integral (see, e.g., ^37] ). 
As follows from Cameron-Martin-Girsanov theorem the 
transformation from Eq. (|B1|I to Eq. (|B5|I can be consid- 
ered as a change of variables ^ ^ C in the path integral 
(El) (see, e.g., H). 

The difference between the solutions (|Bl|l and l(B5|l 
is as follows. The function hj{s,^{t, s)) in Eq. I)B1|I 
explicitly depends on the random velocity field v via 
the Wiener path ^, while the function /ij(s, C(^, s)) in 
Eq. (|B5I) is independent of the velocity v. Trajectories in 
the Feynman-Kac formula ((Bl|l are determined by both, 
a random velocity field and magnetic diffusion. On the 
other hand, trajectories in Eq. (|B5p are determined only 
by magnetic diffusion. Due to the Markovian property 
of the Wiener process the solution l|B5p can be rewritten 
in the form 

h,{t,^) = E{S.,j{t,s,^,X')h,{s,X.')} 

= J Qij{t,s,:xL,x')hj{s,:xL')dx' , (B7) 
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where 

Qij ^1 ^ ) 



[4TTD^{t ~ s)f/^ exp 
X Sij (i, s, X, X ) , 



(x^ - x)^ 

(B8) 



S'y (i, s,x,x') = M^{J{t,s,fi)Gij{t,s,n)} and M/^{-} 
means the path integral taken over the set of trajecto- 
ries /i. which connect points (i,x) and (s,x'). The math- 
ematical expectation E{-} in Eq. (|B7p denotes the av- 
eraging over the set of randorn points X' which have a 
Gaussian statistics (see, e.g., [S^). We used here the 
following property of the averaging over the Wiener pro- 
cess E{Mfj,{-}} — M^{-}. We considered a random ve- 
locity field with a finite renewal time. Due to a re- 
newal, the velocity field as well as its functionals hj{s, x') 
and Qij{t, SjXjX.') in the two time intervals are statisti- 
cally independent. Now we make a change of variables 
(x, x') (x, x' = z-t-x) in Eq. HB7|) . i.e., Qij{t, s, x, x') = 
(3ij(i, s, X, z + x) ~ Qij (t, s, X, z). The Fourier transfor- 
mation in Eq. ljB7|) yields 

hi{t,x.) = J J Qij{t,s,x,k)exp{ik- z)dk 

X J ft,j(s, q) cxp[iq • (z + x)] dqdz . 

Since (5(k + q) = {2tt)^^ J exp[i(k + q) • z)] dz, we obtain 
that 

hi{t,x) = {2TTf J Qij{t,s,yi,-q)hj{s,q) 

X exp(iq • x) dq . (B9) 
In Eq. l|B9p the function Qij{t, s, x, — q) is given by 

Qij{t,s,x,-q) = {2ny^ J Qij{t,s,x,z) 

X exp(iq ■ z) dz . (BIO) 



Substituting (5ij(i, s. x, x') = (5y (t, s,x, z) in Eq. ljB7|l 
and taking into account that x' = z -|- x we obtain 

hi{t,x) = J Qij{t, s,x,z)hj{s,z + k) dz . (Bll) 

Equation ljB10|l can be rewritten in the form 

{2nfQij{t,s,x, - q)cxp(iq-x)= J Qij{t, s,x,z) 

X exp[iq • (z + x)] dz . (B12) 



The right hand sides of Eqs. (|B11|) and (|B12|I coincide 
when h(s, z + x) = e exp[iq • (z -I- x)], where e is a unit 
vector. Thus, a particular solution IjBlip of Eq. (|23l) with 
the initial condition h(s,x') = e exp(iq • x') coincides 
in form with the integral ljB12p . On the other hand, a 
solution of Eq. is given by Eq. ljB5|) . Substituting 



the initial condition h(s, ^) = e exp(iq • ^) = e exp[iq • 
(x + {2D„,y/^w)] into Eq. (ESJ we obtain 

h,{t,x) = M^{J{t,s,C)G^,{t,s,C)eJ 

X cxp[iq • (x + {2D„^)^/^w)]} . (B13) 
Comparing Eqs. (|B11|) - (|B13|) we get 

g,,(i,s,x,-q) = {2^T)-^M^{J{t,sX)G^J{t,s,C) 

xexp[^(2i?„0'/'q•w]} . (B14) 



Now we rewrite Eq. (jB14ll using Feynman-Kac for- 
mula ljBl|) . The result is given by 

g,j(i,s,x,-q) = i2TTy^M^{G,j{t,s,^{t,s)) 

xexp[iq.r]}, (B15) 



where |* = | - x. Substituting Eq. (|BT5|) into Eq. l|B9|) 
we obtain 

h,{t,x) = J M^{Gij{t,s,$,) exp[iq- C]hj{s,q)} 

X exp(iq ■ x) dq . (B16) 

The Fourier transformation in Eq. ljB16|l yields Eq. IjAip . 
The above derivation proves that the assumption (|B2p is 
correct for a Wiener path ^. 



APPENDIX C: IDENTITY USED FOR 
DERIVATION OF EQS. AND SMk 

For the derivation of Eqs. H27|) and (|34|l we used the 
following identity 



iki I /.y(k - -Q, K - Q)f/p(Q) exp(iK-R) dKdQ 



1 
2 

i ( dfij 



l(VsV,C7p) 



dks 



4:\dks 



(CI) 



To derive Eq. IjCljl we multiply the equation 
V • u = [written in k-space for Ui{ki — Q)] by 
Mj(k2)i7p(Q) exp(iK-R), integrate over K and Q, and 
average over ensemble of velocity fluctuations. Here 
ki = k + K/2 and k2 = -k + K/2. This yields 

Jt(^h + ^K, - Q,^ (u,{k + - Q)u,(-k + Ik)) 
xi7p(Q)exp(iK-R)dKdQ==0. (C2) 

Next, we introduce new variables: ki = k + K/2 — Q, 
ica = -k + K/2 and k = (ki - k2)/2 = k - Q/2, K = 
ki + k2 = K — Q. This allows us to rewrite Eq. (|C2|) in 
the form 

^(h + ^K,-Q^k{k-^Q,K^q)Up{Q) 



X exp(iK-R) dKdQ = 



(C3) 



12 



Since IQI < |k| we use the Taylor expansion Vp[/,,(k, R)i7p(R)] - J iKp[f,,{k, R)f7p(R)]K 

/y (k - Q/2, K - Q) ~ (k, K - Q) x exp (iK-R) dK . (C5) 

g,+0(Q2). (C4) 



ia/,,(k,K-Q)^ , ^,^2. 



2 afc. 

We also use the following identities: 



[/,j (k,R)C/p(R)]K = J /y (k,K - Q)C/p(Q)dQ , Therefore, Eqs. yield Eq. iHTl . 
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